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Abstract 


We assessed the structural complexity of vegetation, birds and soundscapes in two areas of Braulio Carrillo Na- 
tional Park, Costa Rica, with different road exposure, as baseline information for conservation and management. 
We measured vegetation structure complexity, surveyed birds and recorded soundscapes in an area adjacent to 
the National Highway 32 (Quebrada) and another area 20 km away (Ceibo). Quebrada presented a more open 
forest structure and lower density of trees and shrubs; lower evenness and higher acoustic complexity, higher 
bioacoustic activity and sound pressure level; higher bird abundance and richness. Ceibo showed a higher density 
of trees and shrubs, higher complexity of vegetation structure, higher proportion of biophonies than anthropho- 
nies, and an acoustic community with higher diversity and entropy. The study sites differed in characteristics of 
their bird communities in an unexpected way; these differences appear to be mostly influenced by the structure 
and complexity of the vegetation surrounding the road. However, the soundscapes seemed to properly reflect 
the habitat condition, and are promising for evaluating the ecological condition of a site. 
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Resumen 


Evaluamos la complejidad estructural de la vegetación, las aves y los paisajes sonoros en dos áreas del Parque 
Nacional Braulio Carrillo, Costa Rica, con diferente exposición a carreteras, como información de base para la 
conservación y el manejo. Medimos la complejidad estructural de la vegetación, registramos las aves y grabamos 
los paisajes sonoros en un área adyacente a la Carretera Nacional 32 (Quebrada) y otra área a 20 km (Ceibo). Que- 
brada presentó una estructura forestal más abierta y menor densidad de árboles y arbustos, menor equitatividad 
y mayor complejidad acústica, mayor actividad bioacústica y nivel de presión sonora y mayor abundancia y ri- 
queza de aves. Ceibo mostró una mayor densidad de árboles y arbustos, mayor complejidad de la estructura de 
la vegetación, mayor proporción de biofonías que antrofonías, y una comunidad acústica con mayor diversidad y 
entropía. Los sitios de estudio variaron en las características de sus comunidades de aves de manera imprevista; 
estas diferencias parecen estar mayormente influenciadas por la estructura y complejidad de la vegetación que 
rodea la carretera. Sin embargo, los paisajes sonoros parecieron reflejar adecuadamente la condición del hábitat 
y son prometedores para evaluar la condición ecológica de un sitio. 


Palabras clave. Aves. Bosque húmedo tropical. Carretera. Ecoacústica. Índices acústicos. Ruido. 
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Introduction 


Roads and their traffic can affect wildlife over large 
areas and, in regions with dense road networks, may 
influence a high proportion of the ecological landscape 
(Cooke et al., 2020a). In Costa Rica, for example, 41.2 % 
of the protected areas are directly or indirectly influen- 
ced by the effect of roads, and 532.52 km of roads are 
adjacent to, or less than 2 km away from protected areas 
(Arévalo & Blau, 2018). 


Due to the complexity of biotic communities in tropical 
ecosystems, roads generate a soaring pressure in their 
vicinity (Irombulak & Frissell, 2000). Among the nega- 
tive impacts related to roads are: a) isolation of wildlife 
populations and road kills (Keller & Largiader, 2003); b) 
proliferation of exotic species along the roadside (For- 
man & Deblinger, 2000); c) changes in the microclima- 
tic conditions at forest edges (Pohlman ef al., 2009); d) 
exposure to chemicals pollutants (dust, heavy metals, 
nutrients, ozone and organic molecules; Laurance et al., 
2009); e) invasions of hunters, miners, colonists and land 
speculators, which can be a cause of additional hunting 
pressure on wildlife (Laurance ef al., 2009); and f) noise 
pollution (Francis & Barber, 2013). 


Many studies have associated road disturbances with 
noise, since it can interfere with the acoustic commu- 
nication, due to a masking effect from vehicular noise 
(Dooling & Popper, 2007). The masking effect occurs 
when the distance at which an acoustic signal can be 
heard is reduced by the anthropogenic background noi- 
se (Marten & Marler, 1977). Many bird species modify 
their vocalization behavior to overcome the masking 
effect by modifying the amplitude, frequency, duration, 
timing or rate of their vocalizations (Pieretti & Farina, 
2013; Roca et al., 2016). Furthermore, noise can impo- 
se selective constraints and pressures on species using 
acoustic signals to communicate and reproduce (Pieretti 
& Farina, 2013). For instance, short-term noise exposu- 
re may decrease nestling body size and increase in the 
oxidative stress of Tree Swallows (Injaian et al., 2018). 


Species spatial distribution can also be affected by road 
exposure. For example, bird movements, abundance, 
occurrence, and richness have been negatively correla- 
ted with road exposure (Reijnen et al., 1997; Laurance 
et al., 2004; Herrera-Montes & Aide, 2011; Cooke et al., 
2020b). In Costa Rica, the number of biophonies (me- 
asured as the seconds occupied by bird vocalizations) 
was significantly lower near the road in Carara Natio- 
nal Park (Arévalo & Blau, 2018). This difference was 
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probably related to decreased bird density where road 
noise was highest (Arévalo & Newhard, 2011). 


Other authors have used acoustic indices to charac- 
terize biotic communities relative to road disturban- 
ces, as they offer a relatively fast method of processing; 
acoustic data (Sueur et al., 2008). For example, Pieretti 
& Farina (2013) used the Acoustic Complexity Index 
(ACI) (Pieretti et al., 2011) to quantify bird vocalization 
dynamics in a landscape influenced by road noise; they 
found a positive relationship between ACI and distan- 
ce to the road. In Costa Rica, Tenez (2016) also found 
higher ACI values at the edge of the road in Santa Rosa 
National Park and attributed this result to bird vocali- 
zations. In addition, sites located 10 m away from the 
road showed the highest proportion of anthrophonies, 
while sites located 500 m away presented the highest 
proportion of biophonies. 


The Braulio Carrillo National Park (BCNP) is one of the 
most extensive protected areas in central Costa Rica, 
with 47 689 hectares. It holds approximately 6000 species 
of plants and thousands of animal species (Boza, 1992). 
This protected area is crossed by National Highway 32, 
through which 80 % of the country’s import and export 
goods are transported (Fernandez, 2013). Furthermore, 
plans exist to expand Highway 32, which might increa- 
se vehicular flow through this protected area. Updated 
information about possible influences of Highway 32 
on biodiversity is required to assess present and future 
impacts, and to support management actions. 


The objective of this study was to assess the struc- 
tural complexity of vegetation, soundscape and bird 
communities in two areas of BCNP with different an- 
thropogenic alteration related to road exposure, as ba- 
seline information for park biodiversity conservation 
and management. The area near the Highway 32 was 
expected to show higher habitat disturbance, higher 
sound pressure level and more disturbed vegetation 
structure than the remote area. In addition, we expec- 
ted also the bird community to be less diverse, the 
soundscape to have a lower proportion of biophonies 
than anthrophonies, and lower acoustic diversity or 
acoustic complexity. 


Materials and methods 


Study area. The BCNP is located in the Central Volca- 
nic Mountain Range in Costa Rica. It receives almost 
6000 mm of precipitation annually, with two brief dry 
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periods between January-April and September-October 
(Blake & Loiselle, 2000). The park belongs to the Phyto- 
geographic Unit of the Caribbean Foothill, characteri- 
zed by abundant and heterogeneous vegetation, and an 
understory dominated by palms (Zamora, 2008). The 
park is crossed by the national Highway 32; this road 
was constructed in 1977 to connect the Central Valley 
with the Caribbean coast. The construction of the road 
implied the opening of forest cover in a radius of up to 
1 km for the passage of machinery and storage of ma- 
terials (Pérez et al., 2009). To avoid further destruction 
of the forest, the national park was created in 1978. 


The study was conducted at two sites within the 
BCNP, each selected based on its distance from 
Highway 32 (Figure 1). The sites are named “Quebra- 
da” (10%09'42.35”N- 83°56'16.11” W; on Highway 32 
and “Ceibo” (10°19'14.25 ”N- 84°04’28.87”W; ~ 20km 
from Highway 32, and ~10 km from another paved 
road). Quebrada is a Ranger Station dedicated mostly 
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to public use, and it is located on Highway 32. The 
area surrounding the trails belonged to cattle farms in 
the mid-20th century, thus the vegetation composition 
still maintains the characteristics of a regenerating fo- 
rest (Pérez et al., 2009). Ceibo is a Sector of the BCNP 
devoted exclusively to protection and research and it 
is located at about 20 km away from Quebrada. The 
access to Ceibo is by a low-traffic unpaved road, and 
there is a minimum distance of 1.5 km from this road 
to the sampling points. Both Ceibo and Quebrada sites 
are included in the life zone “Very Humid Tropical Fo- 
rest Transition to Premontane” and present an altitude 
between 500 and 600 m a.s.l. (TEC, 2014). 


Sampling points: 


At each study site we systematically placed 12 sam- 
pling points in transects along all available trails, where 
permission was granted by park rangers (Figure 1). 
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Figure 1. Location of Ceibo and Quebrada sites within the Braulio Carrillo National Park, Costa Rica, in reference to their dis- 


tance from Highway 32. A, global view of the park; B, zoom view of Ceibo study site; C, zoom view of Quebrada study site 


Figura 1. Ubicación de los sitios Ceibo y Quebrada dentro del Parque Nacional Braulio Carrillo, en referencia a su distancia a la carretera 32. 


A, vista global del parque; B, vista ampliada del sitio de estudio de Ceibo; C, vista ampliada del sitio de estudio de Quebrada 
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We placed the recorders 50 m away from the trail, ona 
central tree selected trying to avoid obstacles that could 
interfere with sound propagation. All sampling points 
were located at a distance of 200 m from each other, to 
reduce the superposition of sounds between two adja- 
cent recording points and ensure a good spatial sam- 
pling of the soundscape (Farina et al., 2011). However, 
we did not quantify detection spaces of the sites (as 
recommended in Darras et al., 2015). 


We made four visits to each study site between June 
2017 and August 2018. Sampling periods were not si- 
multaneous at both sites; however, we systematically 
alternated visits between sites, to get a representative 
sample of the bird and soundscape events that occur 
throughout the year. Two visits were conducted during 
the rainy season of 2017 (Quebrada: June and October; 
Ceibo: September and November), one visit during the 
dry season of 2018 (Ceibo: February; Quebrada: April), 
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and one visit during the rainy season of 2018 (Ceibo: 
May; Quebrada: August). 


Habitat Survey: 


We measured the vegetation structure around each of 
the 12 sampling points. From a central tree we establi- 
shed four 3x20 m plots, to measure the stem diameter 
of trees greater than 5 cm of diameter at breast height 
(DBH). Within each plot, we established two 3x3 m 
subplots to measure the DBH of those individuals be- 
tween 1 and 5 cm with an approximate height 21.30 m; 
otherwise, they were only counted for that diameter 
class. Within these subplots, we established two 2x1 
m subplots to count individuals with less than 1 cm 
of DBH. With the vegetation data, we calculated three 
vegetation structure variables for each sampling point 
(Table 1). 


Table 1. Description of the environmental variables, bird parameters and acoustic indices calculated for each sampling point at Ceibo and 


Quebrada sites of Braulio Carrillo National Park, Costa Rica 


Tabla 1. Descripción de las variables ambientales, parámetros de aves e índices acústicos calculados para cada punto de muestreo en los 


sitios Ceibo y Quebrada del Parque Nacional Braulio Carrillo, Costa Rica. 


Environmental 


Description 


Number of individuals more than 5 cn DBH 
Number of individuals 1-5 cn DBH 


Measures the opening of the forest in relation to the number of opening pixels 
and forest cover for pictures taken at 45° 


Measures the opening of the forest in relation to the number of opening pixels 
and forest cover for pictures taken horizontally. 


Vegetation vertical structural complexity 


Vegetation horizontal structural complexity 


Number of individuals of all bird species. 


Shannon diversity (Shannon & Weaver, 1949). 


Variable Code 
Basal Area basal_total Total Basal Area 
Abundance 5 cm abund_ind_5 
Abundances 1-5 cm abund_1 to 5 
Vertical canopy GAP v 
opening 
Horizontal canopy GAP h 
opening 
Vertical Mean informa- MIG v 
tion gain 
Horizontal Mean infor- MIG h 
mation gain 
Bird parameter Description 
Bird abundance Abundance 
Bird richness Richness Number of species. 
Shannon index Shannon 
Pielou index Pielou 


Acoustic Indices 
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Evenness of Pielou (Pielou, 1984). 


Description 
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Environmental "Ee 
Variable Code Description 
Acoustic Complexity ACI Measures the variability of sound intensity over time and frequency. Higher 
Index values indicate more complex soundscapes (Pieretti et al., 2011). 
Summarizes sound intensity distribution applying the Shannon index over 
Acoustic Diversity ADI the proportion of signals above an intensity threshold across the spectrum. 
Index High values indicate that sound intensity is spread evenly across frequency 
bands (Villanueva-Rivera et al., 2011). 
PE E CET Same as the ADI, but the Gini index is applied across all frequency bins 
index AFI (Villanueva-Rivera et al., 2011). High values indicate sound intensity is res- 
tricted to a narrow frequency. range. 
Is a function of both amplitude and number of occupied frequency bands 
between 2-11 kHz. It is calculated as the area under each curve included all 
Bioacoustic Index BIO frequency bands associated with the dB value that was greater than the mi- 
nimum dB value for each curve. The higher values indicate greater disparity 
between bands (Boelman et al., 2007). 
l ; Shows the median of the amplitude envelope normalized by the maximum 
Median Amplitude : ae 
MAE value to give results between 0-1. Higher values reflect noisier soundscapes 
Envelope 
(Depraetere et al., 2012). 
Calculates the ratio (biophony — anthrophony) / (biophony + anthrophony). 
Normalized Difference NDSI Where biophony is the energy on the frequency range 2-8 kHz and anthro- 
Soundscape Index phony between 1-2 kHz. (Kasten ef al., 2012). Higher values indicate a higher 
proportion of sound intensity in the biophony band. 
Counts the number of major frequency peaks obtained on a mean spectrum 
Number of Peaks NP scaled between 0 and 1. Higher values indicates a greater frequency spectral 
complexity (Gasc et al., 2013). 
Divides intensity values by the sum of intensity in a frequency band. Calcu- 
lates the entropy for these values using the negative sum of each value multi- 
Spectral Entropy Hf plied by the log of the value and divided by the log of the number of intensity 
values (Sueur et al., 2008). Higher values indicate evenly distributed sound 
across frequency bands. 
Divides intensity values by the sum of intensities in a time step. Calculates 
A Ento Ht entropy for these values using the negative sum of each value multiplied by 
P PY the log of the value and divided by the log of the number of intensity values. 
(Sueur et al., 2008). Higher values indicate evenly distributed sound over time. 
It is the product of Hf * Ht. Quantifies the variation across time and frequency 
Total Entropy TE (Sueur ef al., 2008). Higher values indicate evenly distributed sounds across 
frequency and time. 
Mean Sound Pressure SPLMean Calculates the average sound pressure level in dB, using a decibel reference 


Level value of 2x10-5 Pa. Higher values indicate louder soundscapes. 


Additionally, we took four horizontal and four verti- 
cal photographs (at about 45° angle) at each sampling 
point, following the recommendations provided in 
Martin and Proulx (2016), with the modification that 
no photograph touched the ground. For this, we used 
a Canon 80D camera with an 18-135 mm canon lens 
(EF-S f / 3.5-5.6 IS USM, Canon Inc.) and the following 


settings: variable exposure, 15 mm focal length, f / 6.3 
aperture, variable ISO, 2816-1880 resolution, automatic 
focus, image stabilizer on, and JPG type of compres- 
sion. Each group of photographs was sorted by study 
site and sampling point to characterize the geometry of 
the forest using the following variables: MIG (average 
information gain, a measure of structural complexity; 
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Proulx & Parrott, 2008) and GAP (opening fraction, 
which measures the opening of the forest in relation 
to the number of opening pixels and forest cover; 
Rosin, 2001). The estimation parameters were obtai- 
ned by using the “gap_fraction” and “calculateHis- 
to” functions of the “LAI and imagemetrics” package 
(https: / / github.com/cmartin/ImageAnalysisPrimer). 
Four vegetation structure complexity variables were 
obtained for each sampling point (Table 1). 


Bird Community Survey. We conducted bird point 
counts on each sampling point at a distance of about 
20 meters from the recorders, to avoid interference du- 
ring recording times. We conducted two daily coun- 
ting sessions (four counts per visit) during peak hours 
of bird activity (~ 6:00-8:00 and ~ 14:00-16:00). At each 
point, we recorded number of species and number of 
individuals of each species detected visually and au- 
rally during 6 minutes in a 25 meters radius. We ob- 
tained four bird “parameters” for each sampling point 
(Table 1). However, these parameters do not repre- 
sent actual estimated parameters at the population or 
community level, as we did not consider detectability. 
For example, richness represents observed richness 
(number of species) at each sampling point. Abundan- 
ce represents observed total abundance, i.e. the total 
numbers of individuals of all species detected at each 
sampling point. 


Soundscape Survey. We located 12 “Song Meter Di- 
gital Field Recorders 2 Plus” (SM2 +; Wildlife Acous- 
tics Inc.) to make soundscape recordings. We installed 
recorders at the central tree of each sampling point at 
an approximate height of 1.30 m. In each recorder, two 
omnidirectional microphones were placed, so the recor- 
ding was made through two channels (in stereo). Audio 
files were recorded in wav format (16-bit, 44.1 kHz) 
and stored on 64GB capacity SDHC memory cards. 
The recorders were programmed to make continuous 
recordings during bird activity peaks (4:00-6:59 and 
15:00-17:59), and for periods of 10 minutes at the be- 
ginning of each hour between 7:00 and 14:59. This re- 
cording scheme was implemented for two consecutive 
days during each visit to each study site, for a total of 
four separate visits per site. 


The soundscape recordings were previewed aurally and 
visually through a spectrogram to facilitate identifica- 
tion of undesired noise such as heavy rain in the audio 
file, using the Adobe Audition CC v.6 program. We re- 
moved the minute files showing heavy rains, keeping in 
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total 81 060 minutes of recordings used in the analysis 
(42 480 minutes at Ceibo and 38 580 minutes at Quebra- 
da). We analyzed soundscape recordings using eleven 
of the most commonly used acoustic indices, because 
they represent different components of the soundscape 
(Table 1). The acoustic indices ACI, ADI, AEI and BIO 
were calculated using a maximum frequency of 11 000 
Hz. All acoustic indices were calculated for one-minute 
recording fragments and averaged for each study site, 
sampling point and time of day (4:00 to 17:00). 


Data Analysis. To compare environmental variables, 
bird parameters and acoustic indices between sites, 
we conducted Welch tests for each of them, separa- 
tely. This approach was chosen to take into account 
variance heterogeneity of the data. The results were 
represented graphically, showing P value, using an 
alpha of 0.05. 


To conduct a multivariable description of variability at 
both sites, based on the measured variables, we used a 
multiple factor analysis (Pages, 2002). Variables were 
assembled into three groups: environmental variables, 
bird parameters and acoustic indices (Table 1). All va- 
riables, except site, were numerical and therefore scaled 
to have a mean of 0 and a variance of 1 (each variable 
was subtracted from its mean and divided by its stan- 
dard deviation). 


All analyses and graphs were performed using the R 
programming environment (R Core Team, 2019). To cal- 
culate acoustic indices, we used the following packages 
for R: Sinax (https://rdrr.io/ github/ osoramirez/Si- 
nax/), Soundecology (Villanueva-Rivera & Pijanowski, 
2016) and Seewave (Sueur et al., 2016). We used the Bio- 
diversity® package (Kindt, 2015) to analyze bird point 
count data; gestatsplot for comparison of variables be- 
tween sites (Patil, 2018), and FactoMine® for multiple 
factorial analysis (Lê et al., 2008). 


Results 


Comparison of variables between sites: Ceibo showed 
a higher density of trees with a DBH > 5 cm, and shrubs 
with a DBH of 1-5 cm, and higher complexity of forest 
structure (higher MIG, both vertically and horizonta- 
lly). On the contrary, Quebrada showed more open ve- 
getation structure, with larger gaps both vertically and 
horizontally. However, total basal area was similar at 
both sites (Figure 2). 
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Figure 2. Comparison of environmental variables between Ceibo and Quebrada sites in the Braulio Carrillo National Park, 
Costa Rica. Graphic shows Welch test results, 95% confidence interval, p-value (for an alpha value of 0.05) and n. 


Figura 2. Comparación de variables ambientales entre los sitios Ceibo y Quebrada en el Parque Nacional Braulio Carrillo, Costa Rica. El grá- 
fico muestra los resultados de la prueba de Welch, el intervalo de confianza del 95%, el valor p (para un valor alfa de 0.05) y n. 
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Based on bird point counts, we found that bird com- 
munity was more abundant, rich and diverse at Que- 
brada than at Ceibo; however, it was more even at 
Ceibo than at Quebrada (Figure 3). We detected 127 
bird species in Ceibo (33 families and 111 genera), 
and 157 species in Quebrada (40 families and 124 ge- 
nera). Of the total 170 species, 73 (43 %) were shared 
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by both sites, 75 (44 %) were detected only in Que- 
brada, and 22 (13 %) were found only in Ceibo. Bird 
community at Ceibo was mainly composed of forest 
and forest-edge species, while Quebrada had seve- 
ral species from open areas as well (those preferring 
environments with reduced forest cover or early suc- 
cessional states). 
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Figure 3. Comparison of bird parameters between Ceibo and Quebrada sites in the Braulio Carrillo National Park, Costa Rica. 


Graphic shows Welch test results, 95% confidence interval, p-value (for an alpha value of 0.05) and n. 


Figura 3. Comparación de parámetros de aves entre los sitios Ceibo y Quebrada en el Parque Nacional Braulio Carrillo, Costa Rica. El gráfico 


muestra los resultados de la prueba de Welch, el intervalo de confianza del 95%, el valor p (para un valor alfa de 0.05) y n. 


Acoustic indices showed a community with higher 
complexity and bioacoustic activity in Quebrada than 
Ceibo; however, Ceibo presented higher proportion 
of biophonies than anthrophonies, and an acoustic 
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community with higher diversity, higher frequency 
entropy, and total entropy than Quebrada. Indices 
SPLmean and MAE also showed higher sound pres- 
sure level in Quebrada than in Ceibo (Figure 4). 
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Figure 4. Comparison of acoustic indices between Ceibo and Quebrada sites in the Braulio Carrillo National Park, Costa Rica. 


Graphic shows Welch test results, 95% confidence interval, p-value (for an alpha value of 0.05) and n. 


Figura 4. Comparación de índices acústicos entre los sitios Ceibo y Quebrada en el Parque Nacional Braulio Carrillo, Costa Rica. El gráfico 
muestra los resultados de la prueba de Welch, el intervalo de confianza del 95%, el valor p (para un valor alfa de 0.05) y n. 


Multivariable description of variability at the study 
sites: Ceibo and Quebrada sites separated completely 
in the multiple factor analysis; besides, each of the sam- 
pling points clustered around the ellipse that represents 
the corresponding study site (Figure 5). 


The first axis was mainly explained by bird indices 
(richness and abundance) and the vertical opening of 
the forest. The second axis was mostly explained by 
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environmental variables representing abundance of 
shrubs 1-5 cm DBH and abundance of trees > 5 cm of 
DBH (Figure 6). The variables with the most explana- 
tory value are those representing vegetation structure 
and complexity, and bird parameters. However, axis 
one and axis two together explain 47 % of the varian- 
ce, meaning that there are still other factors explai- 
ning data variability, other than the ones measured 
in this study. 
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Figure 5. Multiple factor analysis in Ceibo and Quebrada sites at the Braulio Carrillo National Park of Costa Rica, according 
to 3 groups of variables used (environmental variables, acoustic indices, and bird parameters). Each sampling point is named 
on the graphic using 4 digits, the first two indicate the study site (CE = Ceibo and QG = Quebrada) and the last two refer to 
a consecutive number for each point from 1 to 12. 


Figura 5. Análisis factorial múltiple en los sitios Ceibo y Quebrada del Parque Nacional Braulio Carrillo de Costa Rica, de acuerdo con 3 
grupos de variables utilizadas (variables ambientales, índices acústicos y parámetros de aves). Cada punto de muestreo se nombra en el 
gráfico con 4 dígitos, los dos primeros indican el sitio de estudio (CE = Ceibo y QG = Quebrada) y los dos últimos se refieren a un número 
consecutivo para cada punto del 1 al 12. 
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Figure 6. Multiple factor analysis in Ceibo and Quebrada sites of the Braulio Carrillo National Park of Costa Rica, according 
to three variable groups used (environmental variables, acoustic indices, and bird parameters). The variables with the larger 
value of contribution (darker color of arrow), contribute the most to the definition of the dimensions. Variables that contri- 
bute the most to Dim.1 and Dim.2 are the most important in explaining the variability in the data set. 


Figura 6. Análisis de factores múltiples en los sitios Ceibo y Quebrada del Parque Nacional Braulio Carrillo de Costa Rica, de acuerdo con 
tres grupos variables utilizados (variables ambientales, índices acústicos y parámetros de aves). Las variables con mayor valor de contribu- 
ción (color más oscuro de flecha) son las que más contribuyen a la definición de las dimensiones. Las variables que más contribuyen a Dim.1 
y Dim.2 son las más importantes para explicar la variabilidad en el conjunto de datos. 


BioTA COLOMBIANA 22 (1) - 2021 | 67 


DOI: 10.21068/c2021.v22n01a04 


Discussion 


Our observations suggest that Ceibo was a more ma- 
ture and preserved site than Quebrada, due to its 
higher horizontal and vertical complexity of vegeta- 
tion structure, and higher density of trees and shrubs 
(Oosterhoorn & Kappelle, 2000; Proulx & Parrott, 2009). 
However, Quebrada presented a higher richness, abun- 
dance and diversity of birds than Ceibo, contrary to 
what we expected. This is inconsistent with other stu- 
dies, which have found higher richness and abundance 
in sites with higher complexity of vegetation structu- 
re (Schulze et al., 2004; Mulwa et al., 2012); however, 
other authors have shown a higher bird diversity in 
secondary habitats, which provide a variety of niches 
for species with different habitat requirements (Becker 
& Agreda, 2005; Gray et al., 2007). Furthermore, our 
results are also inconsistent with other studies showing 
a negative correlation between abundance, occurren- 
ce and richness of birds and levels of vehicular noise 
(Reijnen et al., 1997; Herrera-Montes & Aide, 2011; Aré- 
valo & Newhard, 2011). 


We found bird species commonly inhabiting open 
areas, such as Dives dives, Myiozetetes granadensis, and 
Tyrannus melancholicus, only in Quebrada. These re- 
sults could be related to a more disturbed vegetation 
structure in Quebrada, as a possible consequence of 
anthropogenic factors, such as road exposure, tourism 
activity, and historical management in the site (Pérez et 
al., 2009). Roads themselves create open gaps, favoring 
microclimatic edge gradients (Pohlman ef al., 2007), 
and the proliferation of plant species that attract and 
sustain a higher diversity of birds in these areas (Levey, 
1988; Flores & Dezzeo, 2005). Moreover, the nearby 
presence of the Rio Sucio canyon could also contribute 
to explain the higher diversity of birds in Quebrada, 
since it connects the site with grasslands in the lower 
plains of the Caribbean slope (Gillies & Clair, 2008). 


ACI and BIO were higher at Quebrada than at Ceibo, 
like bird parameters. ACI has shown mixed results 
regarding its ability as a bird indicator. Some studies 
found positive relationships between ACI and num- 
ber of bird vocalizations (Pieretti et al., 2011), species 
richness and bird abundance (Eldridge et al., 2018; Re- 
tamosa et al., 2018), or bioacoustic activity in general 
(Towsey et al., 2014). Other studies found weak or ne- 
gative relationships between ACI and species richness 
or biophonic activity (Mammides et al., 2017; Eldridge 
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et al., 2018). Nevertheless, ACI has been more consis- 
tently associated with vehicular noise in other studies, 
showing a positive relationship with traffic noise (Pie- 
rreti & Farina, 2013; Tenez, 2016). These authors ar- 
gued that birds may be attempting to propagate their 
signals with greater emphasis (e.g. amplified redun- 
dancy or song intensity) to offset noise masking effect. 
On the other hand, BIO has been tested in different 
environments, showing positive correlation with bird 
richness (Eldrigde et al., 2018), abundance (Retamosa 
et al., 2018), number of vocalizations (Pieretti et al., 
2011) and biotic diversity (Fairbrass et al., 2017). 


Despite higher bird diversity in Quebrada, this site 
seemed to be less acoustically diverse than Ceibo, with 
sound intensity spread unevenly across frequency 
bands (Villanueva-Rivera et al., 2011, Bradfer-Lawrence 
et al., 2019). Characteristics of some bird species found 
in Quebrada could help to understand these results. 
For example, many species in open areas rely more 
on visual signals compared to species that inhabit the 
interior of the forest, which rely primarily on acoustic 
signals (Laverde-R. et al., 2017, 2018). In fact, Laver- 
de-R. et al. (2017) found evidence that birds that occu- 
py closed or darker habitats have a greater variety of 
song syllables and less colorful plumages. 


Furthermore, we observed bird flocks more frequent- 
ly in Quebrada than in Ceibo, and they also inclu- 
ded more individuals. Some species within the flock 
use special vocalizations to maintain group cohesion 
(Powell, 1985). However, these sounds are generally 
simple with notes designed for close-range communi- 
cation to reduce predator detection risk (Greenberg, 
2001). It has also been hypothesized that mixed flock 
species may have converged using similar vocaliza- 
tions to mimic members of other species (Powell, 1985; 
Goodale & Kotagama, 2006). Based on this, presence of 
more abundant mixed flocks in Quebrada might help 
explain high bird abundance in this site, without ne- 
cessarily resulting in higher acoustic diversity. 


Presence of Highway 32 seemed to reflect in the acous- 
tic environment. For example, SPLMean and MAE in- 
dices indicated an environment with higher sound 
pressure in Quebrada than Ceibo. Moreover, there 
was a higher proportion of anthrophonies than biopho- 
nies in Quebrada than in Ceibo. This is consistent with 
Tenez (2018), who found more anthrophonies near the 
paved entrance road to Santa Rosa National Park, Costa 
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Rica. However, we found positive, high values in the 
NDSI index in both Ceibo and Quebrada, which indi- 
cated that soundscapes at both sites were mostly com- 
posed of biophonies. 


In general, composition of both acoustic communities 
in Ceibo and Quebrada appears to indicate a better 
partitioned acoustic community in Ceibo. According 
to the acoustic niche hypothesis, soundscape is a limi- 
ted resource, and individuals must compete to com- 
municate effectively, resulting in sound niches, where 
the soundscape is spatially and temporally divided 
(Krause, 1987; Pijanowski et al., 2011). Therefore, the 
more complex the habitat, the more complex and spe- 
cific the acoustic partition (Pieretti & Farina, 2013). 


Some authors have criticized the use of population or 
community parameters as indicators of environmental 
disturbance (Temple & Wiens, 1989; Niemi et al., 1997; 
Campos-Cerqueira et al., 2019). For example, abun- 
dance of bird species fluctuates in a specific way, and 
may not represent the trend of other species (Venier 
& Pearce, 2004). Furthermore, changes in bird popu- 
lations are generally associated with various factors, 
making it difficult to establish a cause-effect relations- 
hip with respect to the environmental changes evalua- 
ted (Temple & Wies, 1989). However, further insight 
could be attained by studying particular bird species, 
as well as other parameters related to communication 
behavior, condition, or reproductive success. 


Overall, we could argue for an indirect effect of the road 
on the variables and parameters considered, through 
changes in the structural complexity of the vegetation 
surrounding it. Although bird parameters behaved 
unexpectedly in this study, soundscapes did seem to 
reflect well the habitat condition of both study sites, as 
seen in other studies (Do Nascimento et al., 2020), and 
are promising for evaluating the ecological condition 
of a site. 
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